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ABSTRACT 



Aims. Very few molecular species have been detected in circumstellar disks surrounding young stellar objects. We are carrying out 
an observational study of the chemistry of circumstellar disks surrounding T Tauri and Herbig Ae stars. First results of this study are 
presented in this note. 

Methods. We used the EMIR receivers recently installed at the IRAM 30m telescope to carry a sensitive search for molecular lines in 
the disks surrounding AB Aur, DM Tau, and LkCa 15. 

Results. We detected lines of the molecules HCO + , CN, H 2 CO, SO, CS, and HCN toward AB Aur. In addition, we tentatively detected 
DCO+ and H 2 S lines. The line profiles suggest that the CN, HCN, H 2 CO, CS and SO lines arise in the disk. This makes it the first 
detection of SO in a circumstellar disk. We have unsuccessfully searched for SO toward DM Tau and LkCa 15, and for C-C3H2 toward 
AB Aur, DM Tau, and LkCa 15. Our upper limits show that contrary to all the molecular species observed so far, SO is not as abundant 
in DM Tau as it is in AB Aur. 

Conclusions. Our results demonstrate that the disk associated with AB Aur is rich in molecular species. Our chemical model shows 
that the detection of SO is consistent with that expected from a very young disk where the molecular adsorption onto grains does not 
yet dominate the chemistry. 



Key words. stars:formation-stars: individual (AB Aur, DM Tau, LkCa 15) - 
tary systems: protoplanetary disks 
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1. Introduction 



Circumstellar disks are complex systems in which essentially all 
the processes that play a role in the interstellar medium, UV radi- 
ation, X-rays, grain surface chemistry, molecular depletion, tur- 
bulent mixing, accretion flows and time dependency, are work- 
ing. Chemical models with increasing complexity have been de- 
veloped in the last decade (see e.g. Aikawa et al. 2000; Dutrey et 
al. 2007; Agundez, Cernicharo & Goicoechea, 2008; Nomura et 
al. 2009), but the disk chemistry is a quite unexplored field from 
the observational point of view. Large millimeter telescopes have 
started to provide some insight into the chemistry of the cold gas 
toward the most massive nearby disks. Thus far, few molecules 
(CO, 13 CO, CN, C 2 H, HCN, HNC, HCO + , H 2 CO) have been de- 
tected in circumstellar disks. This small molecular inventory is 
mainly due to the weakness of the molecular emission from cir- 
cumstellar disks. Disks have small masses, lower than 0.1 M , 
small sizes, radii of a few 100 AU, and because of depletion in 
the midplane and/or photodissociation in the surface, the disk 
averaged abundances of most molecules (including CO and its 
isotopologues) are a factor of 5-10 lower than in the interstellar 
medium. High sensitivity is, therefore, required for an observa- 
tional study. We have carried out a sensitive search for molecular 
lines mainly in the disk around the Herbig Ae star AB Aur using 
the IRAM 30m telescope. Some lines have also been searched 
toward DM Tau and LkCa 15. Our results show the rich molec- 
ular content in the disk around AB Aur. 



AB Auriga is one of the nearest, brightest and best studied 
Herbig Ae stars. It has a spectral type A0-A1 (Hernandez et al. 
2004) and is located to the Southwest of the molecular cloud 
L1517 (Duvert et al. 1986), at a distance of 145 pc (van den 
Ancker et al. 1998). Interferometric observations at millimeter 
wavelengths detected the circumstellar disk around this star in 
the continuum and in the CO (and its isotopologues) lines (Pietu 
et al. 2005). Instead of being centrally peaked, the continuum 
emission is dominated by a bright, asymmetric (spiral-like) fea- 
ture at about 140 AU from the central star. 

The disk modeling of the continuum and molecular emission 
showed that the disk is warm and showed no evidence of CO 
depletion. Schreyer et al. (2008) searched for emission of the 
HCO + l-»0, HCN l-»0, CS 2->l, C 2 H l-»0 and some CH 3 OH 
lines in this disk using the Plateau de Bure Interferometer 
(PdBI), but they only detected the HCO + l->0 line. 



2. Observations 

The list of observed lines and the telescope parameters are given 
in Table 1. The observations were done in two observing pe- 
riods, September, 2009 and March, 2010, with the new EMIR 
receivers arranged to provide a bandwidth of 4 GHz in both, the 
3mm and 1mm bands. As backends we used the wide bandwidth 
autocorrelator WILMA which provides a spectral resolution of 
2 MHz and covers the whole band, and the narrow bandwidth 
correlator VESPA centered at the line frequency and providing 
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Table 1 : List of targeted lines 



Line ~ Freq.(GHz) HPBW(") rj b 



CO 2->l 230.538 10 0.63 

HCO + l->0 89.188 28 0.81 

HCO +1 3-^2 267.558 9 0.53 

HCN l->0 88.631 28 0.81 

HCN 3^2 265.886 9 0.53 

CN l->0 113.490 22 0.81 

CN 2->l 226.874 10 0.63 

CS 2->l 97.981 25 0.81 

CS 3^2 146.969 16 0.74 

C 2 H 3^2 262.004 9 0.53 

H 2 CO 3 03 ->2 02 218.222 11 0.63 

SO 3 4 ^2 3 1 38.178 17 0.74 

SO 2 5 6 ^4 5 219.949 11 0.63 

c-C 3 H 2 2 12 ^1 0>1 85.339 29 0.81 

c-C 3 H 2 6 ,6^5i5 217.822 11 0.63 

c-C 3 H 2 6i 6 ^5 05 217.822 11 0.63 

DCO + 2-Ul 144.077 16 0.74 

DCN 2->l 144.828 16 0.74 

SiO 2-»l 86.847 29 0.81 

SiO 6^5 260.518 9 0.53 

HCO l ,i 3/2, 2->0 .o 1/2,1 86.671 29 0.81 

HCO l ,i3/2, l-»0 00 1/2,0 86.708 29 0.81 

HCO lojl/2,1-^0.0 1/2,1 86.777 29 0.81 

HCO l ,i 1/2, 0-»0 .o 1/2,1 86.805 29 0.81 

H 2 S 1i,q->1q,i 168.763 14 0.74 

1 Observed with the wide band spectrometer WILMA in the HCN 3— >2 
tuning. 

2 Observed with the wide band spectrometer WILMA in the c-C 3 H 2 
6o,6->5i i5 tuning. 

a spectral resolution of 80 kHz at 1.3mm and 40 kHz at 2.7mm 
(~0.1 km s _1 ). All the observations were done using the wobbler 
switching (WS) procedure with a throw of 120". This procedure 
provides flat baselines which are essential for detecting weak 
and wide lines toward compact sources, which is the case for 
the lines arising in circumstellar disks. In the case of AB Aur, 
the disk is still immersed in the parent cloud whose emission ex- 
tends farther than the wobbler throw (see Semenov et al. 2005). 
Then, at the velocities of the ambient cloud the detected emis- 
sion is just the ON-OFF balance without any physical interpre- 
tation (remind that the OFF position is moving in the sky during 
the source tracking). For this reason, we have blanked the chan- 
nels corresponding to the ambient cloud emission in the spectra 
toward AB Aur. We have searched for C-C3H2 and SO also to- 
ward DM Tau and LkCal5. In these cases, contamination from 
the ambient cloud is not expected. The observational results are 
shown in Table 2. 

3. Results 

Fig. 1 shows some of the spectra observed toward AB Aur. The 
lines from the molecular cloud are very narrow, Av~0.5 km s _1 , 
and centered at 5.9 km s _1 (Duvert et al. 1986 and Fig. 2a). 
The emission of the ambient cloud lies at the velocities 
[5.4,6.4] km s _1 . The channels corresponding to these velocities 
are blanked in the spectra shown in Fig. 1. 

After blanking the cloud velocities, we have detected emis- 
sion at >3cr of the HCO + l-»0, CN l->0, H 2 CO 3o,3^2 , 2 , SO 
34— >23, CS 3— >2, and HCN 3— >2 lines. In addition, we have ten- 
tatively detected (~3cr) the DCO + 2->l andH 2 S li,o->lo,i lines. 
All the (>3cr) detected lines have the typical two-horn profile ob- 
served in the lines coming from the circumstellar disk, with two 
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Fig. 1. Spectra obtained with the 30m telescope toward AB Aur. 
Dashed lines indicate the velocities at which the disk emission 
arises ([4.2,5.6] km s _1 and [6.5,7.25] km s" 1 ). 
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Fig. 2. a) Spectra of the C ls O l->0 toward AB Aur observed 
with the IRAM 30m telescope by Fuente et al. (2002). We have 
adopted this profile as a pattern profile for the ambient cloud 
emission. Vertical lines indicate the velocity interval [5.6,6.4] 
km s^ 1 . b) Interferometric spectra of the 12 CO 2— »1 toward the 
star position (Fig. 2 (bottom) of Pietu et al. 2005). Note that the 
disk emission occurs at the velocity intervals, [4.2,5.6] km s _1 
and [6.4,7.25] km s _1 . These velocities are indicated by vertical 
dashed lines. 



peaks centered at 4.8+0.25 km s _1 and 6.8+0.25 km s _1 (see 
Fig. 1 and Fig. 2b). This prompts us to interpret the emission 
of these lines as arising from the circumstellar disk. The only 
doubtful case is the CN 1 — >0 line in which the two-horn profile 
is not so clear. Since CN is one of the most abundant species in 
disks (Dutrey et al. 1997, Thi et al. 2004, Oberg et al. 2010), we 
decided to keep it in our list of detected species. The narrowness 
of the CN 1 — >0 line could be due to the fact that its emission is 
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Table 2: Observational results 
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Fig. 3. Comparison between the spectra of the HCO+ 1 — >0 line 
observed toward AB Aur with the 30m telescope and the PdBI. 
Dashed lines indicate the velocities at which the disk emission 
arises ([4.2,5.6] km s" 1 and [6.5,7.25] km s" 1 ). 



coming from the outermost part of the disk. We would like to 
remind, however, that this detection requires of further confir- 
mation by interferometric observations. 

The disk was previously detected in the HCO + 1— >0 line us- 
ing the PdBI by Schreyer et al. (2008). Therefore, we can use 
this line to check the validity of our interpretation. In Fig. 3, we 
compare our HCO + spectrum with that observed toward the star 
position by Schreyer et al. (2008). Since the synthesized beam 
of these observations was 5.2" x 4.8", this spectrum missed the 
emission of the outer part of the disk (R>378 AU). The emission 
of this outer region is expected at velocities <0.87 km s _1 from 
the systemic velocity. We only consider velocities >0.45 km s _1 
relative to the systemic velocity, therefore the outer part of 
the disk is not relevant in our comparison. The integrated in- 
tensity emission of the 30m spectra in the velocity intervals 
[4.2,5.4] km s _1 and [6.4,7.25] km s _1 is lower by a factor of 
« 1.3 than the integrated emission of the HCO + 1— >0 line as ob- 
served with the PdBI (see Fig. 3). The agreement is acceptable 
and we consider that our interpretation is valid. 

In Table 2, we show the list of non-detections (in the consid- 
ered velocity ranges). Our 3<x upper limit to the emission of CS 
2— >1 line improve by a factor of 2 the previous one obtained by 
Schreyer et al. (2008) using the PdBI. 



4. Averaged molecular abundances in the disk 

In the following we derive approximated average column den- 
sities in the disks assuming optically thin emission, a uniform 
temperature of Tj.= 10 K and 20 K, and Local Thermodynamic 
Equilibrium (LTE). The assumed disk sizes (diameters) are: 13" 
for AB Aur (Pietu et al. 2005), 6" for LkCal5 (Thi et al. 2004) 
and 13" for DM Tau (Pietu et al. 2007). In Table 3 we compare 
the obtained fractional abundances with those derived in other 
disks following a similar procedure. The first result is that the 
molecular abundances measured toward AB Aur are very similar 
to those found toward other disks which reinforce our interpre- 
tation of the disk origin for the observed lines. 

The HCO + abundance in AB Aur is similar to those mea- 
sured in the TT stars LkCal5 and TW Hya, and in the HAe stars 
HD 163296 and MWC 480. Only DM Tau presents a signifi- 
cantly larger (a factor of 10) HCO + abundance which suggests 
that DM Tau is a special case among circumstellar disks. The 
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same remains true for CN and HCN. Both molecules present 
abundances of ~10" 10 (CN) and -lO" 11 (HCN) in all the disks 
except DM Tau in which the measured abundances are a fac- 
tor of 10 larger. The case of H2CO is a bit different. It is also 
overabundant in DM Tau (5xl0~ 10 ) but it is underabundant in 
TW Hya (<10~ 12 ). Our SO detection in AB Aur is the first one 
in a circumstellar disk and ours, the first estimate of the SO abun- 
dance in a circumstellar disk. We have unsuccessfully searched 
for SO toward DM Tau and LkCal5. Our upper limit toward 
DM Tau shows that in contrast with the behavior observed in 
the other species, SO is not overabundant in DM Tau relative to 
AB Aur. The abundance of C2H in AB Aur is 1000 times lower 
than in DM Tau. We have used the C2H 3— >2 line to derive the 
upper limit to the C2H abundance. Recent results by Henning et 
al. (2010) show that the excitation temperature of C2H could be 
lower than 10 K in some disks. In this case, the value of the up- 
per limit would increase. The searched C-C3H2 line has not been 
detected in any of the observed disks. The obtained upper limit 
to the [c-C3H2]/[HCO + ] ratio is still consistent with the observa- 
tional results in PDRs (Fuente et al. 2003). 

In Table 3, we also compare some representative column 
density ratios. These ratios are meaningful providing that the 
molecules arise in the same region that could not be the case. 
The value of [CN]/[HCN] is fairly uniform (within a factor of 2) 
among the observed disks. However, there are large variations, 
more than an order of magnitude, in the values of [CS]/[SO], 
[H 2 CO]/[HCO + ] and [C 2 H]/[HCO + ]. This suggests that these 
column density ratios are more sensitive to the in most cases 
poorly known disk structure and/or grain properties. One im- 
portant parameter in disks is the deuteration degree. Qi et al. 
(2008) derived a [DCO + ]/[HCO + ] ratio of -0.03 in TW Hya. 
Our tentative detection of DCO + in AB Aur would imply a 
[DCO + ]/[HCO + ] ratio of -0.03, similar to TW Hya. 
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Table 3: Comparison with other disks 1 
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1 Abundances relative to H 2 have been calculated assuming a disk diameters of 13" for AB Aur (Pietu et al. 2005), 6" for LkCal5 (Thi et al. 2004) 
and 13" for DM Tau (Pietu et al. 2007). The assumed disk-averaged molecular hydrogen column densities are: N(H 2 )=l.lxl0 22 cm -2 for AB Aur, 
2.7xl0 22 cm -2 for DM Tau and 1.4xl0 23 cm" 2 for LkCal5. In our calculations for DM Tau and LkCal5, we adopt T rot =10 K. Abundances for the 
other disks/molecules are taken from Dutrey et al. (1997) (for DM Tau), Thi et al. (2004) and Guilloteau et al. (2006). 



5. Discussion 

In order to guide our interpretation of the observed features and 
provide additional support to their disk origin, we have per- 
formed a (preliminary) chemical model adopting the disk and 
stellar parameters from Schreyer et al. (2008) and the updated 
chemical network of Agundez et al. (2008). Our aim is to in- 
vestigate if detectable SO column densities can be produced in 
this disk. In Fig. 4 we show the radial distribution of the verti- 
cal column densities of some molecules as calculated at 2.5 Myr 
(the age of AB Aur). In agreement with the observational re- 
sults of Pietu et al. (2005), the intense stellar radiation makes 
the disk to be moderately warm, with temperatures above 20 K 
even in the disk midplane so that volatile molecules such as 
CO are not severely depleted on grain surfaces. In our model, 
the major gas phase reservoir of sulphur are the CS and SO 
molecules, with high column densities of SO mainly present in 
the inner R<200 AU region of the disk. The SO abundance de- 
crease rapidly with time because of the adsorption onto the grain 
surfaces. The youthness of the AB Aur disk could be key to have 
higher SO abundance. 

6. Summary and conclusions 

We have taken advantage of the high sensitivity of the EMIR re- 
ceivers recently installed in the IRAM 30m telescope to make 
a sensitive search for molecular emission in three prototypical 
disks, AB Aur, DM Tau and LkCa 15. Our results and conclu- 
sions can be summarized as follows: 

- We have detected the HCO + 1-^0, CN 1^0, H 2 CO 
3o,3->2 ,2, SO 3 4 ^2 3 , CS 3-»2, HCN 1^0 and HCN 3^2 
lines toward AB Aur. In addition, we have tentatively de- 
tected the DCO + 2— >1 and H2S li,o— ^lo.i lines. Based on 
the lines profiles, we interpret the emission of the CN 1 — >0, 
HCN 3^2, H 2 CO 3 ,3^2 ,2, CS 3->2, and the SO 3 4 ^2 3 
lines as arising from the disk. If confirmed, this is the first 
detection of SO in a circumstellar disk. 




100 , v 500 1000 



Fig. 4. Vertical column densities of various molecules as a func- 
tion of radius at 2.5 Myr as calculated by the chemical model. 

- We have unsuccessfully searched for SO toward DM Tau and 
LkCa 15. The obtained upper limits show that SO is under- 
abundant in DM Tau relative to AB Aur. 

- We have searched for C-C3H2 toward AB Aur, DM Tau and 
LkCa 15. The obtained upper limits are still consistent with 
the [c-C 3 H 2 ]/[HCO + ] values obtained in PDRs. 

Our observational work has significantly increased 
(from 1 to 6) the number of species detected toward the 
disk in AB Aur. If confirmed by interferometric observations, 
the SO detection would be the first one in a circumstellar disk. 
Our chemical model suggests that the high SO abundance 
derived in AB Aur disk is consistent with that expected in a very 
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young and warm disk where depletion of gas onto grains is not 
dominating the chemistry yet. 
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